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^ I A review is given of recent theoretical work on the superfiuid dynamics of 

0\ • trapped Bose gases at finite temperatures, where there is a significant frac- 

tion of non- condensate atoms. One can now reach large enough densities and 
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(^ I collision cross-sections needed to probe the collective modes in the collision- 

dominated hydrodynamic region where the gas exhibits characteristic super- 
fluid behavior involving the relative motions of the condensate and non- 
condensate components. The precise analogue of the Landau- Khalatnikov 
'^ ■ two-fluid hydrodynamic equations was recently derived from trapped Bose 

G . gases, starting from a generalized Gross- Pitaevskii equation for the conden- 

Q I sate macroscopic wavefunction and a kinetic equation for the non- condensate 

atoms. 

PACS numbers: 03.75.Fi, 05.30.Jp, 67.40.Bz, 67.40-w. 

1.. INTRODUCTION 

Superfiuid behavior is the most striking property of liquid ^He. We 
recall that spectacular experiments of Kapitza as well as Allen and Misener 
in 1938 first showed that liquid ^He below the transition temperature of 
2.19K could exhibit flow through thin channels without any viscosity. Since 
then, a dominant theme of research on liquid ^He has been to understand the 
origin of this superfluidity and to work out how it modifies the dynamical 
response functions and excitations of liquid ^He. 

By the early 1960's, a successful field theoretic formalism was 
developed,a'EI with the broken-symmetry expectation value of the quantum 
field operator <l>(r, t) = {ip{r, t)) playing the role of the order parameter of the 



A. Griffin and T. Nikuni 

Bose superfluid. In this theory, Bose-condensation was the underlying source 
of all the unique properties of superfluid ^He (equivalence of density fluctua- 
tions with elementary excitations, two-fluid hydro dynamical behavior at low 
frequencies, superfluid flow, quantized vortices, etc). In striking contrast to 
the atomic Bose gases, the explicit role of the Bose condensate is somewhat 
hidden in superfluid ^He. Indeed, this fact has encouraged the development 
of theories for superfluid ^He based on variational groundstate wavefunctions 
(beginning with Feynman in the mid-1950's) without any explicit reference 
to the underlying role of the macroscopic wavefunction describing the con- 
densate. While very useful for computational purposes, these approaches 
have given little insight into the fact that there is a new phase in liquid ^He 
which exhibits superfluidity. For this reason, they do not appear to give 
much hope for developing a general theory of the superfluid dynamics of 
Bose fluids, going from dilute gases to liquid '^He. 

Such a unified description of this kind has taken on a new importance 
as a result of spectacular advances that have occured in BEC research in the 
last year. These discoveries have started to focus research interest in atomic 
gases on issues related to siyperfluidity. The recent creation of quantized 
vortices (at JILA in BoulderH and the ENS in Parisa) have suddenly made 
the BEC community aware of rotating traps, mutual friction, vortex arrays, 
etc, ideas long familiar to the superfluid ^He community. In addition, the 
successful creation of a BEC in ^^Rb gas at JILAB by working close to a 
Feshbach resonance has given us a superfluid Bose gas where the s-wave 
scattering length a can be enormous (o > 10^ A). Already this has moved 
the value of the gas parameter (na^) from the range ~ 10~^ (typical for recent 
BEC experimentsQ using ^^Rb and ^"^Na atoms) to a respectable 10~^. Thus 
atomic condensates are "starting" to overlap on a strongly interacting Bose 
liquid like ^He (where na^ ~ 1, with a now being the hard core diameter). 

More generally, the BEC community is increasingly using the language 
and ideas developed in condensed matter and quantum liquid research. 
Arguments based on topological arguments are used (untwisting of order 
parameters,^ for example). Recent Bragg experiments using second-order 
Raman scattering of light have allowed one to measure^ the dynamic struc- 
ture factor S'(q, u) of trapped Bose gases, a quantity which has played a 
central role in our understanding of superfluid ^He for almost 50 years.0 In 
the last few years, much of the theoretical analysis of atomic condensates has 
come from condensed matter theorists who have worked on superfluid liquids 
(such as Baym, Fetter, Ho, Leggett, Pethick, Pitaevskii and Stringari). We 
hope that the superfluid ^He experimental community will be stimulated by 
the new BEC child in their midst and, in particular, come up with new ideas 
about how to probe the Bose-condensed nature of superfluid '^He. 
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As a contribution to this emerging synthesis, this article reviews recent 
workB'Ej on deriving (from an approximate but still microscopic model) the 
two-fluid hydrodynamic equations for a trapped Bose gas at finite tempera- 
tures. The equivalent two-fiuid equations were first derived phenomenologi- 
cally by Landau in 1941, initially without any reference to a Bose condensate. 
Generalized to include damping from transport coefficients, these coupled 
equations for the superfiuid and normal fiuid components form the basis of 
our understanding of all "superfiuid" behavior in liquid ^He.lliJ This two-fiuid 
description waslater proven to be a consequence of a Bose broken-symmetry 
by BogoliubovEj in 1963, with the superfiuid velocity Vs(r,t) always being 
proportional to the gradient of the phase of <l>(r, i). 

Of course, the two-fiuid hydrodynamic equations describe local equilib- 
rium and require that there be enough collisions to achieve this (cut <^ 1, 
where r is the relaxation time needed to reach local equilibrium). This is 
easy to achieve in a liquid (including liquid ^He) and that is why, histor- 
ically, classical fiuid hydrodynamics was understood a long time before a 
microscopic treatment was available. Using a kinetic equation for interact- 
ing atoms, it was Boltzmann who first (~ 1885) discussed the precise set 
of conditions required for this hydrodynamic description to be valid. Up to 
the present, most experiments on the collective modes of Bose gases have 
probed the opposite "collisionless region" {lot » 1). However one expects 
that the two-fiuid hydrodynamic region will be increasingly studied in the 
next few years, taking advantage of the fact that one can now produce high 
densities of atoms (currently, N ~ 10^ — lO"^ atoms) and also we can enor- 
mously increase the atomic collision cross-section [a = Sna'^) by working 
near a Feshbach resonance (as in recent workQ on ®^Rb). 

As references on BEC research, we recommend the 1998 Varenna Sum- 
mer School Lectured^ as well as the excellent reviewEj by the Trento theory 
group. For a detailed treatment of topics discussed here, see Refs. 14,^. 



2.. DYNAMICS OF A PURE CONDENSATE 

The crucial idea behind the field-theoretic description of Bose con- 
densation developed 40 years agou'B is to isolate the condensate "degree 
of freedom." Thus the quantum field operator is decomposed as '(/'(r) = 
(V'(r)) + ipi^): where $(r) = (V'(r)) describes the Bose condensate (which is 
treated as a classical field). Here V'(i") is the non-condensate component of 
the quantum field operator and satisfies Bose commutation relations. This 
formalism (first developed in a systematic way by Beliaev in 1958) allows 
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one to neatly isolate the condensate. The many-body theory is then devel- 
oped to study the dynamics of the non-condensate fields -0 and -0^ , with the 
condensate playing the role of the "vacuum." A crucial feature is that $(r, t) 
is complex, 



$(r,t) 



nJrAje 



ie(r,t) 



(1) 



with the identification (this is justified later) of the superfluid velocity 
mvs{r) = ?iV6'(r, t). At T = in a dilute Bose gas, one can assume that all 
the atoms are in the condensate $(r,t), which satisfies the famous Gross- 
Pitaevskii time-dependent equation of motionjl^l 
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Here K;x(r) is the parabolic trap potential (a magnetic trap acting on the spin 
of an hyperfine atomic level) and gndv, t) is the Hartree field produced by the 
condensate atoms. At the very low temperatures of interest in BEG studies, 
the atoms have very low energy and one can use the s-wave scattering length 
approximation. In this case, the effective interatomic interaction is given by 
■u(r — r') = g6{r — r'), where g = Airah /m. One should not think of (|2|) 
as simply a Schrodinger equation for a "single-particle" wavefunction, since 
$(r, t) is an order-parameter which is well defined even at finite temperatures 
(see later). As we all have learned over the last 5 years, the T = GP 
equation in (^) contains a huge amount of physics and its solutions have 
dominated BEG research up to nowBllj Many examples are discussed in the 
invited and contributed papers at this QFS Gonference. 

The stationary solutions that the GP equation $o(r, i) = <I>o(r)e~*^''*'^ 
satisfy 
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This can exhibit ground state solutions corresponding to vortices, in 
which there is a time-independent condensate current given by m,VgQ(r) = 
^V^o(r)- In the absence of vortices (irrotational flow), the Thomas-Fermi 
(TF) approximate solution of (^ is very simple and leads to the famous 
parabolic condensate profile 
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This profile is considerably wider than the ground state Gaussian wavefunc- 
tion prediction for a non-interacting trapped gas. 
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It is very convenient to rewrite the GP equation (y) in terms of the 
condensate density nc(r,i) and velocity Vs{r,t) local variables. One finds 

dnc{r,t) / ,N . ,^ 

"i(^ + ^Vv2) = -VMc(r,t), (5) 

where the condensate chemical potential is 

/i^V^ AT 
/ic(r, t) = - ^ ' + Ve^r) + gucir, t). (6) 

Written in this form, the condensate dynamics is "hydrodynamic" looking, 
even in the absence of collisions. This is because the GP equation (|2|) de- 
scribes a large number of atoms in the same single-particle quantum state. 
Linearizing around <l?o(i')) the two coupled equations in (^) can be solved 
for the collective oscillations of the condensate at T = 0. A standard sim- 
plification (first introduced by Stringarill3) is to neglect the "quantum pres- 
sure", the first term in (0). Within this TF approximation, one obtains the 
Stringari "wave equation" for fluctuations in 5nc{r, t) 

^^ = -V . (n,o^ j = ^ V . (n.o(r) V^ne) . (7) 

We note that this wave equation could be equally well written in terms of 
fluctuations of the phase of $(r, t), since 869 /dt = —gSric- As an exam- 
ple of such condensate oscillations, we consider a uniform Bose gas, where 
6nc{r,t) oc g*'*'''""'^*) . This gives the famous Bogoliubov phonon excitations 
of the condensate (first derived in 1947 by a different method) 



iv = cok, Co = W — —. (8) 

V m 

These phonons are physically unrelated to ordinary (hydrodynamic) sound 
waves in a normal fluid. Solutions of (0) corresponding to breathing, dipole, 
quadrupole and surface oscillations of an axially symmetric trap are in ex- 
cellent agreement with experiment (when N > 10^ atoms). Ej'EJ 



3.. TWO-FLUID SUPERFLUID HYDRODYNAMICS 

It is straightforward to extend the preceding T = analysis to finite 
temperatures where there is a large fraction of atoms outside of the con- 
densate described by $(r, t). One then needs two equations of motions (a) 
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A generalized GP equation for <&(r, t) which includes coupling to the non- 
condensate atoms, (b) A Boltzmann kinetic equation for the single-particle 
distribution function /(p, r, t) describing the non-condensate atoms. 

A simple microscopic model has been analyzed in some detail in Ref. ^ 
which appears to contain the correct physics. By resticting oneself to finite 
temperatures, the important non-condensate atoms can be described by the 
simple particle-like spectrum 

ep{r,t) = ^ + Ve^{r) + 2g[nc{r,t)+n{r,t)] 
2m 

This assumes that hioo <^ ksT ; gnc{r,t) <^ ksT, namely both the energy 
level spacing in the harmonic well trap and the average mean field should 
be less than the average kinetic energy of the atoms. Under these condi- 
tions, the high-energy thermal atoms can be described by the single-particle 
distribution function /(p, r, t) obeying a Boltzmann equation 

^I^^ + P .V./(p,r,f)-V.f/(r,t).Vp/(p,r,t) 
ot m 

= C22[/]+Ci2 [/,<!>]. (10) 

This involves two types of collision integrals. Collisions between non- 
condensate atoms are described by 

C22 [/] = T^Its / dP2 / dP3 / dP4<5(p + P2 - P3 - P4) 

X [(l + /)(l + /2)/3/4-//2(l + /3)(l + /4)]. (11) 

In contrast, collisions which transfer atoms between the condensate and non- 
condensate are described by 

Ci2[/, *] = j—y^ j dpi / (ip2 / dp35(mv^ + Pi - P2 - Ps) 

X 5{ec + £pi - £p2 - Sp-i) [5 (p - pi) - 6{p - P2) - 5{p - ps)] 
X [n,(l + /i)/2/3-nJi(l + /2)(l + /3)]. (12) 

Here ec(r, t) = fic + \mv1 is the condensate atom local energy {jic is defined 
below in (|l6|)) and ?tt,v, is the condensate atom momentum. 



Clearly the C12 collisions do not conserve the number of atoms in the 
condensate. Thus they enter (in contrast to C22 collisions) in a direct way 
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in a "generalized" GP equation for <l>(r, t), namel; 



.,5$ 

'""m 



■— + Fex(r) + gnc{r, t) + 2gfi{r, t) - ihR{r, t) 

2m 



$ . (13) 



There is a new term {2gh) associated with the Hartree-Fock mean field of 
the thermal cloud on the condensate. In addition, there is a dissipative term 
associated with C12 collisions, 

^ ' ' J (27r)3 2nc{r,t) 2nc(r,t) ^^ ^ ^ ' 

As with the T = GP equation of Section ^, one can rewrite (13) in terms 
of the nc(r, t) and Vs(r, t) variables to find (compare with (^ and (|^) 

f/T) 

-£ + VncV, = -ri2[/,$] 



V2 /^ 

M,(r, t) = - ^ ^ " + Kx(r) + (7'^c(r, t) + 25n(r, t). (16) 



where now 



We note that even at finite T, the condensate equations of motion in ( |T5|) 
are "hydrodynamic looking" except that now in the continuity equation for 
the condensate, there is a source term ri2 associated with C12 collisions. 

One has to solve these coupled equations for $(r, t) and /(p, r, t) self- 
consistently - and clearly Cu will play a special role. While (y) and ([l3|) give 
a sound basis for discussing the general non-equilibrium behavior of trapped 
gases at finite T, the collision-dominated hydrodynamic region is especially 
interesting since then the non-condensate can be described in terms of a 
few "coarse-grained" local variables. If the C22 collisions are rapid enough 
(relative to the frequency cu of the collective mode of interest, u;t22 ^ 1), 
then /(p,r, t) will be driven to the local equilibrium Bose distribution 

f{p,r,t) = — -^ . (17) 

This function uniquely satisfies C22 [/,*!'] = 0. With this local equilibrium 
distribution, the thermal atoms are now completely described in terms of 
the /3, v„,/i and U variables, all dependent on (r, t). Inserting / given by 
(|l7| ) into the Boltzmann equation (y) and taking moments with respect to 
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the momentum p, one easily derives a set of hydrodynamic equations for the 
non-condensate variables n, v„ and jl. This procedure was first developed 
by Boltzmann for deriving the hydrodynamic equations for a classical gas. 

To summarize, at this stage we have arrived at a closed set of cou- 
pled equationg§ll2Hl3 for the condensate variables {ric'Vs, fJ-c) and the non- 
condensate variables (n,v„,/i). Linearizing these equations around static 
equilibrium, one can derive the hydrodynamic collective modes which corre- 
spond to coherent motions of both the condensate and non-condensate. Our 
new two- fluid equations exhibit a new relaxation time r^, which describes 
how fast the two components come into diffusive equilibrium with each other 
(i.e., how fast fidiffi^jt) = /i(r,t) — fj.c{r,t) relaxes to zero). Of course, in 
static equilibrium, we have flo = ^co but when we perturb the trapped Bose 
gas from equilibrium, one can have 6fiiiiff{T^,t) ^ and it relaxes to zero in 
a time r^. As one would expect, r^ is proportional to the collision time T12 
associated with C12 collisions. 

How is all this related to Landau's famous two-fluid hydrodynamic equa- 
tions? These equations are expressed using different local variables. The 
linearized form of these two-fluid equations areO 



d6 



n 



dt 



-V-Jj 



d5j^ d6P „ dVe, , d 



27? (d^, - ^TrD6^, 



dt dx^ dxfj, dxi, 

d 

+ ^ — [CiV- (mnco(5vs -(5v„) -FC2V- (5v„] 

or 

fn—^ = -V [(5/U + mCsV • mnco(5vs - (5v„) + mC,A^ ■ Svn] 

d6s 1 

— = -V.{so6vn) + -V-{KV6T). (18) 

Here P{r,t) is the local pressure, s(r,t) is the local entropy density (en- 
tirely associated with the normal fluid) , and Dfj,,y = I ( "3J^ + ^^ ) • The 
expressions 

6p = m5n = 5ps + dpn 
mSj = psoSvs + PnoSvn- (19) 

make the two-fluid nature of the theory clear. We have included hydrody- 
namic damping from various transport processes (the thermal conductivity 
K, the shear viscosity r/ and the four second viscosity coefficients Ci, all depen- 
dent on position through the local condensate density nco(r)). In making 
the detailed comparison between the Landau two-fluid equations and our 
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two-fluid equations, a key role is played by the condensate relaxation time 
Tfj_ introduced above. 

In the limit ujt^ — > 0, one can showcl that our two-fluid equations are pre- 
cisely equivalent to the Landau equations summarized in (|T8|), but without 
the dissipative terms. This makes sense, since these equations are derivedEJ 
under the assumption that the superfluid and normal fluid are always in 
local equilibrium with each other. This corresponds to /u(r, t) = ^ci^,t), ie, 
T^ -^ 0. However, even in this limit, it turns out that ri2 in ( [iq ) is finite and 
is crucial to ensure the equivalence with the Landau equations. One finds 
one can make the identification (valid within our model) 

Ps{r,t) = mnc{r,t) 

Pn{r,t) = mh{r,t). (20) 

In the limit of wr^ <C 1,ujt22 *C 1, our equations can be shownEJ to 
reduce precisely to the Landau-Khalatnikov two-fluid equationsEJ in (|18|). 
What is very satisfactory is that, in our model, the second viscosity damping 
terms turn out to be proportional to Sfidiff, i-e., to the fact that /Xc / A- This 
is physically quite reasonable and expected. Thus the four second viscosities 
are all proportional to the collision time ri2 associated with the C12 collision 
term. In contrast, the other transport coefficients {k and 77) are associated 
with deviations for /(p, r, t) from the local equilibrium distribution ( |T7| ) and 
have contributions from both the C22 and C12 collision integrals. 

Finally, one can consider the new limit lot^ ^ l,u;r22 ^ 1, which can 
arisety at temperatures close to Tbec (where r^ becomes very large due 
to a sort of "critical slowing down"). This is outside the region of validity 
of the usual Landau two- fluid equations. Surprisingly, however, our new 
linearized two-fluid equations lead againll3 to the Landau-Khalatnikov two- 
fluid equationsEJ but now with frequency-dependent second viscosities 

o(^) = r^—- (21) 

This form could have been predicted, since it is typical of a situation where 
a fluid is coupled into an "internal" degree of freedom with a long relaxation 
time - in the present case, this is the Bose condensate. 

As one might expect, our new form of the two-fluid equation leads nat- 
urally to a central zero-frequency relaxational mode which is missed by the 
usual LK two-fluid hydrodynamic equations. One can showE-3'tJ that as 
T -^ Tbec from below, this mode (which is strongly coupled to the ther- 
mal conductivity) becomes the usual zero frequency thermal diffusion mode. 
Observation of this mode below Tbec is a goal for future BEC studies since 
this mode is a unique feature of the hydrodynamics of a trapped Bose gas. 
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In conclusion, we have sketched recent progress in deriving the equiv- 
alent of the dissipative Landau-Khalatnikov two-fluid hydrodynamic equa- 
tions for a trapped Bose-condensed gas. The nice feature about our deriva- 
tion is that our microscopic model allows us to compute all the thermody- 
namic functions, transport coefficents and relaxation times which enter into 
these equations. What is needed now is a sustained experimental effort to 
check the predictions of this two-fluid theory. 

ACKNOWLEDGMENTS 

Much of the work presented in Section 3 has been done in close collabo- 
ration with Eugene Zaremba. A.G. is supported by NSERC of Canada and 
T.N. is supported by the JSPS. 

REFERENCES 

1. P.C. Hohenberg and P.C. Martin, Ann. Phys. (N. Y.) 34, 291 (1965). 

2. For a review, see A. GrifRn, Excitations in a Bose-condensed Liquid (Cambridge, 
N.Y., 1993). 

3. K.W. Madison et al, Phys. Rev. Lett. 84, 806 (2000). 

4. M.R. Matthews et al, Phys. Rev. Lett. 83, 2498 (1999). 

5. S. Cornish et al, Phys. Rev. Lett. 85, 1785 (2000). 

6. Bose-Einstein Condensation in Atomic Gases, Proceedings of the International 
School of Physics "Enrico Fermi", ed. M. Inguscio, S. Stringari and C. E. Wie- 
man (LOS Press, Amsterdam, 1999). 

7. M.R. Matthews et al, Phys. Rev. Lett. 83, 3358 (1999). 

8. D.M. Stamper-Kurn et al, Phys. Rev. Lett. 83, 2876 (1999). 

9. E. Zaremba, T. Nikuni and A. Griffin, J. Low Temp. Phys. 116, 227 (1999). 

10. T. Nikuni, E. Zaremba and A. Griffin, Phys. Rev. Lett. 83, 10 (1999). 

11. I.M. Khalatnikov, An Introduction to the Theory of Superfluidity {W. A. Ben- 
jamin, N. ¥., 1965). 

12. N.N. Bogoliubov, Lectures on Quantum Statistics, Vol. 2 Gordon and Breach, 
N.Y., (1970), p.148. 

13. F. Dalfovo, S. Giorgini, L.P. Pitacvskii and S. Stringari, Rev. Mod. Phys. 71, 
463 (1999). 

14. A. Griffin, in 13'^ Physics Summer School: Bose-Einstein Condensation, ed 



by C. Savage and M. Das (World Scientific, Singapore, 2000), in press; cond- 



mat/0006382 



15. S. Stringari, Phys. Rev. Lett. 77, 2360 (1996). 

16. T. Nikuni, private communication; T. Nikuni and A. Griffi n, to be published. 

17. T. Nikuni, A. Griffin and E. Zaremba, |cond-mat/9911336| . 



